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Effective stochastic resonance (SR) is numerically and analytically studied using a model with 
coupled two particles exposed to heterogeneous, i.e., particles dependent, amplitude of noise. Com- 
pared to previous SR models of single particle and to those of coupled two particles exposed to 
equivalent amplitude of noise, the present model exhibits a more intensive resonance of, at least, 
one particle exposed to the non-larger amplitude of noise with the assistance of another particle. In 
a certain range of conditions, this effective resonance of one particle overwhelms the poor resonance 
of the other particle, meaning that heterogeneous amplitude of noise leads the system, not only 
locally but also in the average of the whole, to the effective SR. 



I. INTRODUCTION 



As reviewed irP, various kinds of phenomena have been recognized as stochastic resonance phenomena (SR) and have 
been studies with a variety of methods^^. Among theoretical models for SR, the simplest is that described by the 
Langevin dynamics of a particle confined in a double-well potential periodically deformed by an oscillating external 
field, 

where, 

V{x) = -ax 2 + bx 4 - Acos(Qt)x, (2) 

the last term of which indicates the oscillating external field. Hereafter, we call the above model as the original SR 
model. 

To quantify the efficiency of SR, the signal-to-noise ratio (SNR), the Spectral Power Amplification (SPA) and other 
quantities have been proposed. According to recent relating studies^, here, we employ SPA as the index of SR. SPA 
is, roughly, the power ratio between the input and the output signals, and is defined as, 

J^g S out {uj)duj 

where f2, Si n (u>) and S ou t(oj) represent the frequency of the input signal, power spectrums of input and output 
signals, respectively. The denominator means the power of the input signal (sinusoidal wave), and the numerator is 
the corresponding power of the output signal. Here, S is taken to be larger than the peak width of the power spectrum. 
In our simulation, considering the resolution of the numerical Fourier transformation, ^ is approximated by 

Sin (£2) 

A series of efforts to obtain the higher SPA or SNR have been conducted through variations of the original SR model, 
one of which was made by coupling two particles applied in a common double- well potential periodically deformed by 
the external field as described by Through t he rea sonable tuning of coupling strength, this coupled particles model 
gives a higher SNR than the original SR modeP2EH. As another directional variation, cascade dynamics using one- 
way coupling between oscillating sub-systems was simulated^, that exhibits a high synchronization to the oscillating 
external field by adding a proper amplitude of noise even below Hopf bifurcation threshold. Note that most of the 
previous models of SR have treated systems in which constituent elements are exposed to equivalent amplitude of 
noise, whereas a few models have considered the heterogeneous, i.e., element dependent, amplitudes of noise. As 
discussed in one of such models^, the noise emerging in a living cell is not always equivalent through the cell because 
the origin of noise in living systems is, in addition to the thermal fluctuation, the local density fluctuation of reacting 
molecules in each cell. 

Therefore, as one way for revising the preceding SR models considering at the time the potential comparison 
to actual systems, we study SR model of coupled two elements exposed to the heterogeneous amplitude of noise, 
restricting the model as simple as possible to focus only on the basic mechanism of the underlying system. 
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FIG. 1: Contour lines of ASPA in D 1 - D 2 space for (a) if = 1.00, SI = it /A, (b)K = 5.00, SI = it/ 4, (c)K = 1.00, SI = tt/64, 
(d)K = 5.00, SI = 7r/64, (e)extended figure of the left-lower part of (d). The brighter tone means the higher value of ASPA. 
ASPA is symmetric with respect to the diagonal line Di = D2 from the definition of ASPA. Only in the case of (d)( and (e)), 
the peaks of ASPA fall off the diagonal, equivalent amplitude of noise, line. In (d), contours of ASPA > 0.230 are drawn by 
bolder lines and those of ASPA < 0.222 are abbreviated to focus on the maximum points 



II. LANGEVIN DYNAMICS MODEL AND SIMULATIONS 



The present model is described as, 

dx\ dV(xi) 
dt dxi 

dx 2 _ dV(x 2 ) 
dt 8x2 



+ K(x 2 -x 1 )+£ 1 (t), 

(5) 

+ K(x 1 -x 2 )+&(t), 



where V(x) is same as a — 8.0, b — 0.25, A — 10.0 and K is coupling constant. The quantities £i(t) and £, 2 (t) are 
Gaussian white noise with the relations, 

(&(*)) =0, (6(*)&(s))=2A<M(*-s) (ye{i,2}). (6) 

It should be notified that noise amplitudes Di and D 2 are independently varied as a set of control parameters. 
We numerically calculate ^ and measure SPA for respective particles with varying a set of parameters (Di,D 2 , K). 
The numerical integration is carried out using the 4th order Runge-Kutta method for deterministic part and the 
Euler-Maruyama scheme for stochastic part with time mesh At = 0.005. All data are taken as the ensemble averages 
of 100 simulations with each run performed from t — to t = 102 x 2tt /SI under the randomly chosen initial conditions 
of Xi (xi <G {— yja/2b, -\/a/26}), where we take away the dynamics from £ = 0to< = 2x27r/fHo avoid the dependence 
to initial conditions. In the following, we show the numerical outcomes. Specifically, SPA of particle 1 named SPAl, 
and the average SPA of two particles named ASPA=(SPAl+SPA2)/2, are focused on. 

Fig|l] shows ASPA in D\ — D 2 parameter space for different sets of (K, SI). Note that Fig[lja)-(d) are symmetric 
with respect to the diagonal line D\ — D 2 from the above definition of ASPA. As the oscillation frequency of external 
field, two values, SI = ir/4 and SI = tt/64:, are chosen, and as the coupling strength, K = 1.00 and K = 5.00 are 
chosen. FigjTJa) and (b) show ASPA numerically obtained with SI = tt/4 in the cases of weak (K = 1.00) and strong 
(K = 5.00) couplings, respectively. For both cases, the peaks of ASPA roughly falls on the line D\ = D 2 > meaning 
that equivalent amplitude of noise leads the system to the most effective SR in Di — D 2 space in these cases of 
SI = 7r/4. With a slower oscillation (SI — tt/QA) of external field, if coupling strength is weaker than a critical value, 
the equivalent amplitude of noise still gives the largest ASPA as seen in Figfljc). On the other hand, as shown in 
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FIG. 2: Maximum ASPA obtained in each D1 — D2 space with different K value (o), Maximum ASPA obtained under equivalent 
noise (Di = D2) (°), and corresponding set (D\,D2) (□ and x) for (a) Q = n/4, (b) Q, — 7r/64, here, only the sets satisfying 
Di < D2 are drawn. 



FigjTJd) (and its enlarged view FigfTJe)) for a stronger coupling [K = 5.00) with Q = ir/64, the peaks of ASPA 
in D\ — D2 space are definitely located off the equivalent amplitude line of noise. This means that heterogeneous 
amplitude of noise leads the system to the most effective SR in D\ — D2 space in this case. 

In Fig{2] maximum ASPAs obtained at each if-fixed D\ — D2 space (hereafter max-ASPA) are connected, varying 
K, as graphs. The symbol o in Fig{2ja) shows the relation between K and max-ASPA in the case of O — it /A. Here, 
SR of the system is reinforced as the coupling between particles is strengthened until max-ASPA reaches a peak value 
at finite K (K ~ 2). In the same figure, the symbols, □ and x respectively give the sets (Di, D2) realizing max-ASPA 
at each K. As mentioned above, the value of ASPA is symmetric with respect to D\ and D2 from the definition of 
ASPA, therefore two sets of (Di, D 2 ) realizing max-ASPA are seen if £>i 7^ D 2 - Fig|2l (and Figj4]) show only the sets 
satisfying D\ > Z?2- Two symbols □ and x roughly collapse on each other in FigMa) indicating that max-ASPA 
is attained under equivalent amplitude of noise in this case. Fig{2jh) shows the case of slower oscillation fl = 7r/64. 
Like in Figpja), the symbol shows the relation between K and max-ASPA in each if- fixed D\ — D2 space. In this 
figure, max-ASPA graph consists of two parts divided by a dip around K ~ 3, each of which has a mound shape, and 
as mentioned afterward, the right-side mound including the highest point in the graph is mainly contributed by the 
particle exposed to the lower amplitude of noise. The symbols, □ and x, in the same figure give the set (D\, D2) of 
noise amplitude with which max-ASPA is realized at each K. Obviously different from the case of O = 7r/4 shown 
in Figj2^a) , D\ and D2 in Figj2](b) separate from each other for 2.25 < K < 9.25, within which interval the highest 
value of the max-ASPA is realized. Namely, the most effective SR in D\ — D2 space is realized when heterogeneous 
amplitude of noise is applied to the system. To close up further the effect of heterogeneous amplitude of noise, we plot 
in Fig. 2(b) maximum ASPA obtained along Di = D 2 line in each K-fixed Di — D 2 space. Compared with max-ASPA 
obtained without the constraint D\ = D2 in the same figure, we see the right-side mound in max-ASPA graph is 
definitely ascribed to the heterogeneity of adding noise. 

In the next, we concentrate on the SR of particle 1. Fig{3]ja)-(d) show SPAl (SPA of particle 1) in D\ — D2 space 
for different sets of (K, fi) . The peaks in all of these figures are located off the diagonal line Di = D 2 and within 
the area of D\ < D2. This means that the particle under a smaller amplitude of noise dominantly experiences the 
effective SR caused by the heterogeneous amplitude of noise (hereafter, heterogeneous SR). Figga) and (b) show the 
relation between K and maximum SPA of particle 1 (hereafter max-SPAl) in the cases of = 7r/4 and — tt/64 
respectively, where the sets (Di,D 2 ) realizing max-SPAl in each K fixed D\ — D 2 space also are drawn. In the same 
figures, for the easiness of comparison, fitting curves of the max- ASPA shown in Figj^a) and (b) are added as dashed 
lines, where the graphs of max-SPAl start from the values same as max-ASPA at K — 0. In the small K region, as 
K increases, the values of max-SPAl monotonically increase, and after passing over shoulder parts around which the 
increasing rates once decreases, they recover rapid increasing rate until reaching maximum values. Through whole 
range of K, max-SPAl is kept higher than (or at least same as) the max-ASPA of dashed lines. This, again, means 
that the particles under a smaller amplitude of noise dominantly experiences the effect of the heterogeneous SR. 

If we pay attention to the K value (K > 5.75) realizing the peak of max-SPAl in FigQa) and the decreasing curve 
of max-ASPA around the corresponding K value, the SR of particle 1 is considered not sufficiently significant to be 
reflected in the average SR of the whole system in this case of £1 = 7r/4. Meanwhile, as shown in FigQb), with the 
slower oscillation of 17 = 7r/64, the K value giving the peak of max-SPAl realizes the peak of max-ASPA at the same 
time. Namely, the resonance effect of particle 1 is strong enough to overcome the poor resonance of particle 2 in this 
case of slow oscillation (fi = 7r/64) of external field. In such way, heterogeneous amplitude of noise leads the system 
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FIG. 3: Contour lines of SPA1 in D 1 - D 2 space for (a)K = 1.00, tl = tt/4, (b)K = 5.00, Q = tt/4, (c)K = 1.00, fi = tt/64, 
(d)K = 5.00, = 7r/64. The brighter tone means the higher value of SPA1. 
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FIG. 4: Maximum SPA1 obtained in each K-fixed D\ — D2 space (o), fitting curve of Maximum ASPA (dashed line), and 
corresponding sets {D\,D2) (□ and x) for (a)f2 = tt/4, (b)Q = ir/64. Here, only the sets satisfying Di < D2 are drawn. 



in a certain range of conditions, not only locally but also in the average of the whole, to an effective SR. 



III. ANALYSIS USING MASTER EQUATION MODEL 

To understand the mechanism of the present effective SR of coupled two particles under heterogeneous amplitude of 
noise, we analytically treat the simplified transient dynamics among finite number of states using a master equation; 

dP(a 1 ,a 2 ,t) , . 

dt [ ' 

= W(-a 1 ,a 2 ,t)P(-a 1 ,a 2 ,t) - W(cn, a 2 ,t)P(a 1 , a 2 , t) + W{~a 2 , a 1 ,t)P{a 1 , -a 2 ,t) - W(a 2 ,a 1 ,t)P(a 1 ,a 2 ,t), 

where <7j G {s, — s} represents the state in which particle i € {1, 2} is at the left minima (crj = —s) or the right minima 
(<Xj = s) of the potential expressed by (2) (s = yja/2b). W((Xi, <7j,t) denotes the transition rate from state (c^, <jj) to 
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(— <7j, (Tj), the specific form of which is assumed as, 

W{a h aj, t) = ~ {1 - cx^- tanh(K/A)} {1 - ff< tanh(A cos fit/ DA} (i ^ j), (8) 

ZTi 

where A, K, Di are the amplitude of the input signal, the coupling strength and the noise strength added to the 
particle i, respectively. The quantity is the Arrhenius-type relaxation time, Tj = r exp(A/D i ), where r = 
2tt/ {\ — 2a + /c|(4a + k) and A is the activation energy. This transition rate has a similar form as that introduced in 
the kinetic Ising model for D\ = D^. Here, we should note that this master equation model does not take into account 
the direct transition between states (ci, a 2 ) and (— ax, —02). It indicates that this model has a direct correspondence 
to (5) only at the weak coupling between two particles. 

As shown in Appendix, master equation (7) is analytically solved with several approximations introduced in ref(9) 
for Di = D 2 though some details of calculation are different for the present case of D X ^D 2 . Figj5);a)-(d) show ASPA 
obtained using the analytical solution of (7) in D\ — D 2 space with different sets of (K, £1). 

In Fig[5ja) and (b), the peak of ASPA falls on the line D\ = D 2 > which means the equivalent noise leads the 
system to the most effective SR in D\ — D 2 space. This figures qualitatively corresponds to Figjlja) and (b). With 
slower oscillation O = 7r/64 and weak coupling K = 1.00 the peak of ASPA is still on the diagonal line D\ = D 2 > 
as shown in Fig{5]jc), however the peaks of ASPA fall off the diagonal line with stronger coupling K = 5.00 (Fig{5jd)) 
meaning that heterogeneous amplitude of noise drives the system to the most effective SR in D\ — D 2 space. In 
this way, as long as looking the contours of ASPA, qualitative features of Figjlja)- (d) obtained by the simulation of 
Langevin dynamics model (5) seems to be reproduced by the analytical solution of master equation (7). However 
some remarkable aspects of the heterogeneous SR are not sufficiently described by the present form of master equation 
model as discussed in the followings. 

To discuss the quantitative aspects, max-ASPAs obtained at each if-fixed Di — D 2 space are connected in Figj6] 
as graphs with varying K. Solid lines in Figj6^a), (b) show the relation between K and max- ASPA in the cases 
of il = 7r/4 and f2 = 7r/64, respectively. Here, SR of the system is reinforced as the coupling between particles is 
strengthened until max-ASPAs reach peak values at finite K (K ~ 0.6 in Figj6^a) and K ~ 0.7 in Fig|6jb)), and 
thereafter monotonically decrease. In the same figures, dashed and broken lines respectively give Di and D 2 which 
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FIG. 5: Contour lines of ASPA in D 1 - D 2 space for (a)A' = 1.00, fl = tt/4, (b)K = 5.00, Q = tt/4, (c)A = 1.00, Q = tt/64, 
(d)A = 5.00, Q = tv/64, A = 10.0, To = 2ir/y/(\ — 2a + k\(4a + k), A = 64.0 obtained from the analytical calculation of master 
equation (7). The brighter tone means the higher value of ASPA. ASPA is symmetric with respect to D\ and D2 from the 
definition of ASPA. In the case of (d), the ASPA has two peaks and the peaks fall off the diagonal line. 
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FIG. 6: Maximum ASPA obtained in each K-fixed D\ 
line, -D2: broken line) for (a)f2 = 7r/4, (b)Sl = 7r/64 
r = 2tt/V(| - 2a + k\{Aa + k), A 



D2 space (solid line), and corresponding sets (D\,D2) (Di : dashed 
Here, only the sets satisfying D\ < D2 are drawn and A = 10.0, 
64.0. The curve of D2 in (a) comes out of the graph over a critical value of K. 



collapse on each other in a range of small K and separate beyond f2-dependent critical values of K. The effect of 
the heterogeneous amplitude of noise exceeds that of equivalent noise as K increases. These features shown in Fig. 6 
especially for K < 3 reflect the qualitative characters obtained by Langevin dynamics model (Fig. 2) for the same K 
range. However, it should be cared in Fig. 6(b), that the effect of heterogeneous amplitude of noise is not sharp enough 
to make the second mound of the max-ASPA graph shown in Fig. 2(b) of Langevin dynamics model, nor the overlap 
of Di and D 2 in Fig. 2(b) giving max-ASPA for large K is recognized in Fig. 6(b). Moreover the dominant region of 
heterogeneous SR for large K in Fig. 6(a) has no counterpart in Fig. 2(a). 



IV. DISCUSSION AND SUMMARY 



These distinctions between behaviors of the Langevin dynamics model and the master equation model would 
partially originate from the assumption of weak coupling limit as mentioned in the previous section. However the 
more intrinsic origin for this distinction seems linked to the basic mechanism of the present heterogeneous SR as 
explained below; 

If we rewrite the Langevin dynamics of (5) like, 

dxi _ dV e s{x 1 ,x 2 ) 

dt dx\ 
8x2 _ dVcs{xi.x 2 ) 

dt dx 2 

an overall effective potential V e(xi,x 2 ) for particles 1 and 2 is, 



-&(*). 

(9) 

■&(*), 



V c g(x 1 ,x 2 ) = a{x\ + x 2 2 ) + b{x\ + x\) + K(n - x 2 ) 2 - Acos(£lt)(xi + x 2 ). (10) 

If we take x 2 as constant in the above form, an effective potential of particle 1 becomes monostable above a critical 
coupling strength K c that depends on Qt and x 2 . Then the signal and the coupling force act in the same direction if 
fit = 0, x 2 = s(= \/a/2b), (or if Qt = tt and x 2 = — s(= —y/a/2b)), supposing coupling strength satisfies K > K c , in 
which case the signal is superthreshold. This means that particle 1 will move quickly into the single minimum (if it is 
not already there) without surmounting a potential barrier. Noise is not necessary for this transition of particle 1. On 
the other hand, particle 2 should be sufficiently noisy to trigger the monostable potential for particle 1. In particular, 
the mean escape time of particle 2 should be much smaller than half the period of the signal. In this case, particle 
1 can exhibit essentially a classical (i.e., deterministic) resonance, which results in a large ASPA at unequal noise 
amplitudes. In the present simulation of Langevin dynamics model, the critical coupling strength K c is estimated as 
K c ~ 2.4 using the specific values of parameters: a = 8.0, b — 0.25, A = 10.0. 

This scenario explains the split of max-ASPA graphs over K ~ 2.4 ~ K c in Fig. 2(b), one is max-ASPA within 
whole Z?i — D 2 space and the other is that on Di — D 2 line. Also K c is close to the initiating point of the right-side 
mound of max-ASPA in the same figure. Moreover, it should be noted that, in the limit K — > 00, max-ASPA in 
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Di — D 2 space degenerates on a line D\ + D 2 = constant, causing the vanishment of heterogeneous SR in this limit, 
which is reflected in Fig. 2(b) as the rejoining of D\ and D 2 at large K accompanied with the rejoining of max-ASPA 
values in whole Di — D 2 space and along Di = D 2 line in the same figure. 

To summarize, we numerically showed that the efficiency of SR is enhanced by the coupling of bistable elements 
exposed to heterogeneous amplitude of noise. This heterogeneous SR is caused through a 'task allocation' between 
two particles such that one particle (particle 2) under large fluctuation easily jumps over the potential barrier and 
pulls another particle (particle 1) to the potential minimum state keeping itself under a large fluctuation, then only 
the latter (particle 1) directly contributes to the heterogeneous SR. Therefore if we focus only on the SR of particle 1, 
it is found that the present mechanism works for a wide range in K (coupling strength) and f2 (frequency of external 
signal) as seen in Fig. 4, where the maximum value of SPA1, an SR index of particle 1, is marked under heterogeneous 
amplitude of noise in almost all range of K in the graph. 

However, the most noticeable in the present study is that the effect of heterogeneous SR is largely pronounced 
when coupling strength between two particles exceeds a critical value and the external signal for particle 1 reaches 
super-threshold level with the help of highly fluctuating particle 2. In that case not only SPA1; an index of SR of 
particle 1, but also ASPA; an index of averages degree of SR of particles 1 and 2, mark the maximum value at each- if 
under heterogeneous amplitude of noise. 

In this way, the present heterogeneous SR is an outcome of an exquisite combination of the triggering particle 
with a large fluctuation and the following particle which faithfully responds to the external signal under small (or 
zero) amplitude of noise. We suppose that the heterogeneous SR has potential applicability to a wider variety of 
stochastic resonance phenomena like those of devising artificial sensors with high susceptibility^] anc l that it has a 
certain connection to the basic study on the function of non-isothermal multi-elements systems^. 
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Appendix 



In this appendix, we analytically derive the time evolution of the expectation values of <x; treated in (7). For 
simplification we assume <Ji G {1, —1}. From equation (7), we obtain the dynamics for the expectation values, 



(01O2) = — 
Tl 



tanh (/3ij4cos(fit)) + 771(02) — 771 tanh (/3iAcos(£lt)) (01(72)) 
tanh (p 2 A cos(£lt)) + 772(01) — 772 tanh (p 2 Acos(Slt)) (0102), 

(-(0102) +771 + tanh(/3iAcos(fM))(0 2 ) -771 taah(/3ij4cos(f2t))(<7i)) 



(A.l) 
(A.2) 

(A.3) 



+ — (-(0102) +772 +tanh(/?2Acos(f2t))(oi) - 772 tanh(/3 2 ^cos(r2/j))(0 2 )) . 

T2 



where 77^ = tanh [K Pi) and Pi = 1/-D;. To solve the equation, we assume that input signal is sufficiently weak 
(PiA <C 1), then, the correlation function (ai(t)a 2 (t)} has even parity with respect to A. Consequently, only the 
zeroth order term in A remains in this correlation in the linear approximation. An additional assumption is that 
(ai(t)a 2 (t)} is in a steady state, which means that (ai(t)a 2 (t)) is expressed like, 



(01027 = 



Then (Al) and (A2) are simplified like, 



(02) 



T2 



i_ m 

Tl Tl 



771 T 2 + 77 2 Tl 

n + T 2 



^1; 

<*2> 



Hi cos(fM) 
H 2 cos(fM) 



(A.4) 



(A.5) 



Hi 

Hi 



(1 -77i(0i0 2 ))/3iyl 

3 

Tl 

(1 -772(0102)) P2 A 

T 2 
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The long time limit solution of the above equation is 

(o»(t)) longtime = Rn cos(Qi - tp{) - R i2 cos(Qi - (fi 2 ) = R% cos(Qt - ipi) (ie{l,2}), (A. 6) 

where, 



Ai = * (- (pi + pi) + \J (pi - pi? + ^mmpip^j , 



L = 2?7ipi, 



M 



(Pi - Pi) + - Pif + ^ViViPiPi, 



N = {pi - P2) - \l {pi - P2) 2 + ^mmpipi, 



NHi - LH 2 
(N - M)y/X[+T^' 
-MHi + LH 2 

M(NH 1 - LH 2 ) 



Rn 
Ri 2 
R 2 i 
R 22 



Ri = \/i?-i+^2+ 2i ?il^2COs(^i-^ 2 ), 

ifi = tan — 



ipi — tan 
Pi 



L{N - M)^X( + n 2 

N(-MHi + LH 2 ) 
L(AT-M) v /Al + 2 ' 



l ( Rn sin ipi + R i2 sin ip 2 
Ru cos v?i + R% 2 cos (^ 2 7 ' 



1 

7 

r, = exp(A;i - Xit ). 



Hence, the auto-correlation function is 

fait + At) <Tj(t)) l on gtj me = i? 2 COs(fit + OAt - Vi) COs(fit - ^i), ( A -7) 

and averaging it over initial phase, we get 



((<Ji(t + M)<Ti(t))) = — I ((Tiit + A^aiityion^dt 

Rf cos(ftAt) 



(A.8) 



To obtain the power spectra, we apply the Wiener-Khinchin theorem to (A.8). The power spectrum Si(ui) defined 

by 

/oo 
e-^da^t + At)ai(t)))dAt, (A.9) 
-00 



is calculated in the form 



Si(u) = ^[j(fi-w) + 5(n + u) . (A.IO) 
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And the power spectrum about input signal (Acos(f2f)) is written as 

/oo 
(Acos(nt + nAt)Acos(nt)) exp(-iujAt)dAt 
-oo 

nA 2 

= —(S{n-uj) + s{n + uj)). (A.n) 

By substituting Si for S out in (4), we obtain the local SPA, i.e. SPA,; (i e {1, 2}), and the average SPA, i.e. ASPA, 

SPAi = J, (A.12) 
ASPA = ^ + SPA 2 = ^+3. (A . 13) 

In the present case of a 6 {s, — s}, to comparing with Langevin Equation (5) Ri(K, A) is written as Ri(s 2 K,sA), 
the expectation value cr^ and SPA also are simply rewritten as 

(0i)longtime = sRi(s 2 K, sA) cos(ilt - vjA, (A. 14) 

spa, = ^KM? (A15) 



* Electronic address: nishimor@hiroshima-u.ac.jp 

1 L. Gammaitoni, P. Hanngi, P. Jung, F. Marchesoni, Reviews of Modern Physics, 70 (1998) 223. 

2 R. Benzi, A. Sutera, A. Vulpiani, J. Phys. A, 14 (1981) L453. 

3 J.K. Douglass, L. Wilkens, E. Pantazelou, F. Moss, Nature, 365 (1993) 337. 

4 J. P. Levin, J.P. Miller, Nature, 380 (1996) 165. 

5 J.J. Collins, T.T. Imhoff, P. Grigg, Journal of Neurophysiology, 76 (1996) 642. 

6 E. Simonotto, M. Riani, C. Seife, M. Roverts, J. Twitty, F. Moss, Rhys. Rev. Lett, 78 (1997) 1186. 

7 I. Hidaka, D. Nozaki, Y. Yamamoto, Phys. Rev. Lett, 85, (2000) 3740. 

8 Y. Morishita, T.J.Kobayashi, K.Aihara, J. Theor. Biol, 325 (2005) 241. 

9 R. Takahashi, M. Suzuki, Physica A, 353 (2005) 85. 

Y. Jiang, H. Xin, Phys. Rev. E, 62 (2000) 1846. 

1 A. Neiman, L. Schimansky-Geier, Phys. Lett. A, 197 (1995) 379. 

2 J. F. Lindner, B. K. Meadows, W. L. Ditto, M. E. Inchiosa, A. R. Bulsara, Phys. Rev. Lett., 75, (1995) 3. 

3 P. Jung, P. Hanggi, Phys. Rev. A,44 (1991) 8032. 

4 R. Kawai, A. Awazu, N. Nishimori, preparing. 



